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**F ILE**1D**MTHDMOD 
Lala! MM TTTTTTTTTT «HH HH DDDDDDDD all ‘ala 000000 DODDDDDD 
‘ala MM TTTTTTTTTT HH HH DDDDDDDD MM ala 000000 DDDDDDDD 
MANN 6MMMM TT HH HH DD DD LI 00 00 DD DD 
MMMM =6MMMM TT HH HH DD DD MMMM MMMM 00 00 OD 
a 6M OM TT HH HH DD DD MM MM MM «(OOO 00 DD 
se a TT HH HH DD dD MM 00 00 OD 
Lala! ala! TT HHHHHHHHHH =D DD MM mm 600 00 D0D 
Lala! mM TT HHHHHHHHHH =D DD MM MM 600 00 OD 
Lala! mM TT HH HH DD DD MM MM 600 00 DD 
Lala! MM TT HH HH DD DD MM MM 600 00 DD 
lala! Lala! TT HH HH DD DD MM mm 600 00 DD 
MM Lala! TT HH HH DD DD MM mM 600 00 OD DD 
lala! MM TT HH HH DDDDDDDD MM MM 000000 DDDDDDDD 
mn Lala) TT HH HH DDDDDDDD MM MM 000000 DDDDDDDD 
LL 111111 SSSSSSSS 
LL IIIT] SSSSSSSS 
LL I] SS 
LL I] SS 
LL I] SS 
LL I] SS 
LL I] SSSSSS 
LL I] SSSSSS 
LL I] SS 
LL I] SS 
LL I] $$ 
LL I] $s 
LLLLLLLLLL II111 SSSSSSSS 
LLLLLLLLLL HI! SSSSSSSS 
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Breer obe Tissiee LATANTS Beets lel ns, | hee Cl, 


TOILE BINSBRCD 
IDENT /3-001/ ; File: MTHDMOD.MAR Edit: JCW3001 


—MmPARBARAAAAAAAAAAASAAAAAAALALELALE LAL E SESE SEALERS EEE EE EE ESSER REESE E EEE EEE ES SY 


COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


-* ® 
:® ® 
ad ® 
:® + 
® it 
it THIS SOFTWARE 1S FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED 
# ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE * 
* INCLUSION GF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER 
;* COPIES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY * 
© OTHER PERSON. NO TITLE TO AND OUNERSHIP OF THE SOFTWARE TS. HERESY 
:* TRANSFERRED. . 
** ® 
* THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE * 
;# AND SHOULD NOT RE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT ®* 
:® CORPORATION. . 
® ® 
*® & 
*® * 
:® os 
*® ® 
® * 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 


: LARRARAAAAAAASLAAALALALASASLELALLALELALELELERESAE SERRE ERAS ERASER ESR SESE SESS SS 
: “FACILITY: MATH LIBRARY 

; ABSTRACT: 

This module contains the routine MTHSDMOD: 

It returns the pone inner of the division of argi/arg2 using 


the yen toutes Bg ned 
argl int (argl/arg2)) *arg2 


SOSOSSOSOSSOSOOSOOSSOSOSSOSOSOSOSOSSOSSOSOOSOOSOOSOSOOOOOOSOOS 
SSOSOOSOSSOOSOOSOSOSOSOSOSOSOSOSOSOSOSOSOSOSOSOSOSOSOSOSOSOSOOSOOSOOSOOOO 
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; AUTHOR: Jeffrey C. Wiener, CREATION DATE: 21-DEC-1982 
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NMOS OONANE WIN O OD NAME WN OOD NAME WIN OC OONOU SW OOONOUE Wh 
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00 : 

09 ; MODIFIED BY: 

0 

00 an 

Be -SBTTL HISTORY ; Detailed Current Edit History 
4 ! 3-001 Original version of complete re-write JCW 21-DEC-82 


vr secs ors 
~SBTTL DECLARATIONS 
INCLUDE FILES: 
NONE 


EXTERNAL SYMBOLS: 


GBL 
“EXTEN MTHSSSIGNAL 
“EXTRN MTHSK_FLOUNDMAT 
“EXTRN MTHSK~INVARGMAT 
: LIBRARY MACROS CALLS: 
SSFDEF 


: EQUATED SYMBOLS: 


oo 


“x00 
“XFF 


"Oo 
"oO 
o= 
“Dw 
“oo 
"Oo 


EXP_55 
HIGA_MASK 
OWN STORAGE: 
NONE 
PSECT DECLARATIONS: 
' .PSECT _MTHSCODE 
: CONSTANTS: 


TWO_EXP_55: 
-LONG “x00005C00, “x0 


Sooooooooooo 


NAME WMO ODNAU EWN OOD NOAU EW OO OONOAUSE 


009090909 09 0909 SIN NI NII NN NNN AA AAA AAAOA IIT 


SO 


< 38¢ 91: 3: 7 AX/VMS Macro V04-00 Page (3) 


1: MTHRTL.SRCJMTHDMOD .MAR; 1 


3; Force all external symbols to be declared 


; Define SF$ (stack frame) symbols 


3 55*2°7 


PIC, SHR, LONG, EXE, NOWRT 


3 20055 


—————$ —— + 


v04-00 P | 
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3-001 MTHSDMOD = D REAL*8 remainder 6-SEP=19 1:22: MTHRTL.SRCJMTHDMOD.MAR; 1 


3 - ~SBTTL MTHSOMOD - D REAL*8 remainder | MTH 
4 3 FUNCTIONAL DESCRIPTION: | 
38 3 Return the remainder of argt/args in ats point format 
3¢ 3 Remainder = argl - (int(argl/arg2))*arg | PSE 
8 38 ; The algorithm used to evaluate the DMOD function is as follows: =ee 
0 98 : X = the first argument. MT 
00 3 3 Y = the second ef ounent. 
38 100 ; step 1. m = the exponent of Y. 
101 ; n = the exponent of X. 
B88 1 ¢ s c=n-™m 
0 1 ; If c < 0, end with result = X. 
B38 : step 2. I = the fractional part of Xx. Pha 
43 105 ; J = the fractional part of Y. ove 
000 196 : If lor J, Lele J Ini 
$008 198 ; tep 3 ” 8 toe i, wh 56 for D_floati b Pas 
3 ste ° = 2“(p-l)*I, where p = or oating numbers. 
0008 109: step A. 7s b/s ° Syn 
0608 110 ; T = [T#2*(p-1)J-2*(p-1). T is int(L/J) or int(L/J)+1 Pas 
0008 111; TetLeJjret Sym 
0008 ig 3 ami¢cé,teged T was int(L/J)+1 Pse 
0008 113; step 5. c =c = (p=1) Cro 
0008 114; If ¢ > 0 go to step 3. Ass 
0008 115; step 6. If c = Pe go to step 9. 
0008 116; step 7. L = 2*(p-ltc) * I The 
0008 117; step 8. l=t-J*T 141 
0008 118; step 9. Result = 2*m * | The 
0008 119; 137 
0008 120; 0p 
B008 13 3 CALLING SEQUENCE: 
443 ’ 3 Remainder.wd.v = MTHSDMOD (dividend.rd.r, divisor.rd.r) 
44 138 ; INPUT PARAMETERS: 
000 126 ; Mac 
0008 127; The two input parameters are double precision floating-point ona 
008 128 ; values. They are passed by reference. -82 
00 129 ; 
00 130 0G 

00000004 000 131 DIVIDEND = 4 : Dividend = X in the algorithm. 

00000008 4 : § DIVISOR = 8 : Divisor = Y in the algorithm. The 
008 134 | MAC 
oS 135 ; IMPLICIT INPUTS: 

0 1 $ ; 
00 1 3 NONE 
0 138 3 
173 3 FUNCTION VALUE: 
0 141 ; Remainder of the division of argi/arg2 is returned as a 
3 146 3 double precision floating point value. 
4 144 ; IMPLICIT OUTPUTS: 
00 145 ; 
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oo 
ww 
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7D 
CA 


63 
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snro— 


85 


om } } 00000000 00 00 00 00 0D CD CD CD CD OD 0D 
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SOCSSCSOSSOSOOSOSOOSOSOOSOOSOOOSOSOOOSOSSSSOSOSOSOSOSOOOOSOOOSOOOOOOOOOOOoO 
MRO OONOA NEADS “OC ODNAUN EWN $$ O ODNAUE WIN 0 OONOA UE WW 0 OONOAU EW OOOne 


Ronny @ tt 
DOOOOO OOO OO O00 00960 09 69 09 09 09 09 C9 SI NINN SN NNN NOAA AAA AAA AMI UTUTUIUTI  £ 


WONVHLMW 7 WD WOOCOCOCOCCOCOCCOCOCOCODAA $$ AOP PONS LOLS 


PIP & & EE 


NONE 


NONE 


SIDE EFFECTS: 


COMPLETION CODES: 


AX/VMS Macro 
MTHRTL.SRCIMT 


DMOD.MAR; 1 


MTHS_INVARGMAT = Invalid argument to math Library if the divisor » zero. 


MTHS_ 
th 


PUSHL 


BICw2 
XORW 


BICw2 
XORW 


In STEP_4 and STEP_8 the calculation of 


FLOUNDMAT = Ag beg 
FU bit is set in the callers PSL. 


MTHSDMOD, “M<R2, R3, R4, RS, R6, R7, RB 
@DIVISOR(AP), R2 ; R2/R3 = Y 

; Y=0. Division by zero 
Q@DIVIDEND(AP), RO 3; RO/R1 = X 
#*XFFFFBO7F, R2, R6 ; R6=m is the biased exponent of Y 
#*XFFFFBO7F, RO, RB 3; RB=n is the biased exponent of X 
R6, RB R4 = ¢ = n-m unbiased 
STEP_2 esult is X if X<Y, ie, if c<0 

RO/R1 = X 

R6 ; push m onto the stack 
#*XFFBO, R2 ; R2/R3 = J = unbiased ifract(Y): 
#*x4000, R2 ; J = properly biased ifract(Y)! 
#*XFFBO, RO ; RO/R1 = I = unbiased ifract(Xx): 
#*x4000, RO ; 1 = properly biased ‘fract(Xx): 


underflow in math Library is signaled i 


I=L - J*int(L/J) must be 


computed as precisely as possible. To do this we will need to write J as 
= 

where J1 = the high 24 bits of J and J2 = J - J1, the low 24 bits of J. 

Been OEE is used to extract the 8 bits of J from longword2 that belong 

o JT. 


R2, -(SP) 
#HIGH_MASK, 4(SP) 
(SP), R2, -(SP) 


@DiVIDEND(AP) 


replaces the value 
uh an, the SP 


coce 
mo 


cu 


If I<J 
go to STEP_5 
el = I- 


se Il = 
go to sTeP_5 if 1>0, or 
else the algorithm ends 
the sign of the result is 


v04-00 p | 
M age 4) | 


| 
i 


16-SEP-1984 01:19:04 VAX/VMS Macro V04-00 P 
MTHSDMOD = D al remainder 6-SEP-1 1382 te :22: 50 MTHRTL.SRCJMTHDMOD.MAR; 1 set 3). 
1 BGEQ DONE 3; the same as the sign of | 
50 95 7 MNEGD RO, RO ; the first sroument. A. 


3 DONE: RET 


50 = 01 OF 79 ERROR: ASHQ #15, ; Y=0. Reserved operand 
7E OO'BF 9A MOVZBL antisk TENRARGMAT, =-(SP) t error code 

00000000'GF 01 fe CALLS , G*ATHSSSIGNAL ; signal the error 

50 00001880 8F CO STEP_3: ADDL2 #EXP_55, RO 3 RO/R1 = L = 2e*#(p-1)*] 
3% 
: STEP4; 
; 2*(p=1) = 2*(55) is added and then subtracted from 
: T = int(L/J) to ensure that T = chopped(L/J) or chopped(L/J) +1 


DIVD3 R2, RO ; R6/R7 = T 


L/J 
T+2**(p-1) 


56 550 
56 80 AF 


67 » R6 
60 ADDD2 TWO_EXP_55, R6 3; R6/R7 T 
56 FF7B CF 62 SUBD2 TWO_EXP_55. R6 3; T-2**(p-1) = L/J chopped or choppe 


The calculation of I = L = J*tint(L/J) must be <eees as preciseiy 
as possible. To ,%9 this we will need to write T a 


T=72 
where 71 = the high 24 bits of T and 22 = T = 21, the low 24 bits of T. 


SOoQCOQoooooooQoooQoQooQooooooooo°ooo 
pe eee feololofejojeoleololeoleoleoleoleleololololololeolololo) 
0909008 09 09 0008 SI SISSIES SSNS NOOO AO 
POAAAAA $F OOOOOOOOO ONIN — F SN MO} 


: 
: 
; 
: 
0086 : 
; 
5 
; 


PUPVIVPVPVSUSISIITT BB BS BS ES EE EE ANNAN ANNI AIP NPNINININID 2 9  - S S - OOOOOOO 
WOO NA NE WWN  O OD NAME WIN OC OD NAME WIN $ OC ODNOA UE WWD (OOO NOUS WR "OOONOUW 


0086 : 
0086 ; Now, using J = J1 + J2, 
0086 ; 
0086 : L=-J * int(L/J) = Lb = (1 + J2) © (21 + 22) 
0086 ; = b= (21 * J1) = (21 ® J2) 
0086 : - (22 * J1) _ (22 ® 32) 
0086 : =L-(21*5)2> (224) 
0086 ; 
0086 3° 
0086 
54 56 DO 0086 MOVL R6, R4 : 
55 57 FFFFOFFF 8F CB 0089 BICL3 #H{IGH MASK, R7, RS > RG/RS = 71 
56 654) «662 «(0091 SUBD2—R4, + R6/R7 = 22 
7E O8 AE 54 65 0094 MULD3 R4, B(SP), -(SP) : Compute 21*J1 
20 BE 62 0099 SUBD2 (SP)+, RO > RO/RI = L - gist 
4 6€ 64 009C MULD2 (SP), R4 > R4/RS = 21*%J 
50 54 62 O009F SUBD2 R4, RO : RO/RI = L = 218J 
54 56 OB AE 6 OOA MULD3 8(SP), R6, R4 > RG/RS = 22*J1 
50 54 62 OA SUBD2 = R4, RO : RO/RI = L = 21") = 22*J1 
56 6E) 0 64SCAA MULD2 (SP), R6 > R6O/R7 = 12eJ2 
50 56 62 OOAD SUBD2 -R6 Rg > RO/RI = L - 2*J 
05 14 0080 BGTR STEP : 
6D 13 00B2 BEQL RETURN : End if RO/R1=0 
50 52 60 00B4 ADDD R2, R ; Add J back in because you had 
bat $ T=chopped(L/J)+1 
58 00001880 &F @ STEP_5: SUBL2 #EXP_55, RB rc = co(p-1) = ¢-55 
4 BGEQ P 


+7 AX/VMS Macro V04- 


00 
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58 00001880 8F C0 pce 9 ADDL2 #EXP_55, RB 3c = (p-1)4c = 554¢ 
0C7 $6 3+ 
0C7 65; 
Boer ee : The next two Lines of code are STEP_6 and STEP_7. 
00C7 96 3° 
00C7 6 
so $8 «tb gore | Se8 Rept? RBC RO PL = 1e2*(c4t) 
sb = I*2*(c+ 
OOCC 70 
pect 71 3+ 
Occ ie 3 2*(p-1) = 2*(55) is added and then subtra + i from 
Bore 6 : T = int(L/J) to ensure that T = chopped(L or chopped(L/J)+1 
OOCC 75 
56 50 52 67 OOCC 76 DIVD3 R2, RO, R6 3 ne ls =T=tL/J 
56 «6 FF2C CF «6060: 0000 37 ADDD2 TWO_EXP_55, R6 : R6/ =T= _Teeee (pet) 
56 FF27 CF 62 43 $78 SUBD2 TWO_EXP_55, R6 $ Foote (p}) L/J chopped or choppe | 
OODA 80 ;+ 
OODA 81; 
OODA 36 : STEP_8: 
OODA 283; 
OODA 284 ; The calculation of I = L = J*int(L/J) must be computed as precisel 
Bree +H ; as possible. “Ye, d0 this we will need to write T as 
a +44 3 where Z1 = the high 24 bits of T and 22 = T = 21, the low 24 bits of T. 
OODA 289 ; Now, using J = J1 + J2, 
OODA $39 3 
OODA 91; L=-J* int(L/J) = Lb = (J1 + J2) * (21 + 22) 
OODA $36 3 =L = (21 * J1) = (21 * J2) 
OODA 29 ; - (22 * J1) _ (22 * J2) 
OODA 94; =t-(212eJ) 2 (22 2 J ) 
OODA 295; 
OODA 296 ;:- 
OODA 97 
54 56 DO OQOODA 98 MOVL R6, R4 : 
55 57 FFFFOFFF 8F CB OODD 299 BICL3 M#HIGH_MASK, R7, RS > R4/RS = 21 
56 54 62 OOES 00 SUBD : R6/R7 = 22 
7E 08 AE 54 65 O0E8 01 MULD R4. 8(SP), -(SP) ; Compute Z1*J1 
50 BE 62 ODED 08 SUBD2 (SP)+, RO : RO/R1 =L - gist 
54 ef 64 4 8 MULD Hg », R4 : R4/RS = Z1*J 
50 4 6 F 4 SUBD 0 : RO/R1 = L = Z1*J 
54 56 O08 AE 65 O0F6 305 MULD RCgp) R6, R4 > RG/RS = z2sJ1 
50 54 62 OOFB 06 SUBD2 RG RO. > RO/RI =L- aa ie - 22*J1 
56 6E 064 COOFE 6 MULD2 (SP) > R6/R7 = 
50 56 62 0101 8 SUBD2 = R6 _ : RO/R1 = te - 1*J 
05 14 0104 309 BGTR STEP 9 : 
19 13 1 $ 10 BEQL ig TURN : End if RO/R1=0 
50 52 60 Q1 11 ADDD 3; Add J back in because you had 
198 \ : T=chopped(L/J)+1 
10 AE 4098 8F ae 108 14 STEP_9: SUBW #18000. 16(SP) ; Remove bias from m and 
50 10AE A 111 15 ADDW2 16(S RO : form RO/R1 = 2*meL 
0B 19 «+0115 16 BLSS UNDERF LOW : Branch if underflow 


| 


MTHSDMOD 
3-001 MTHSDMOD = D REAL*8 remainder 
117 17 TEST_SIGN: 
04 f ° 117 + TSTW 
11A 1 BGEQ 
50 8000 8F A 11¢ 0 BISw2 
04 : 1 1 RETURN: RET 
1 ; : UNDERF LOW: 
50 7C¢ 1 4 CLRQ 
0D 04 AD 06 €E1 4 3 2 BBC 
00000000'8F DD 0129 27 PUSHL 
00000000'GF 01 FB big 8 CALLS 
04 0136 9 NO_FU: RET 
tHe 0 
0137 31 END 


D 10 
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@DIVIDEND(AP) : the sign of the result is 
; the same as the sign of 
#*x8000, RO ; the first argument, X. 


RO Set up default result to 0.0 
#SFSV_FU, SFSW_SAVE_PSW(FP), NO_ Fu 
; Branch if caller ase 5 enabled F | 


page FLOUNDMAT : Report MTH$_FLOUNDMA 
, G*ATHSS$SIGNAL : si nal the Weadtaten 
; Return 


MTHSDMOD 6-S 
Symbol table 6-S 


DIVIDEND = 00 4 
DIVISOR = 8 
DONE 


EXP SS = PP 1880 


P=} 1384 1:19: G4 AX/VMS Macro V04-00 
P=1984 11:22: MTHRTL.S RIM MTHDMOD .MAR; 1 


o 
_ 
bob 2) 
oo 
rr 


HIGA_MASK = FFFFOFFF 
RIUWSSS IGNAL etereree xX 0 
THSDMOD 00000008 RG 0 
RTHSK FLOUNDMAT ‘teeeeene x 8 
MTHSK— _INVARGMAT tteeenee x 0 
NO_FU- 0000136 R 8s 
RETURN 4 1R 0 
FSV_ = 00000006 
SFSWISAVE_PSW = 00000004 
STEP_ 0000002A R 8s 
STEP_ 000000 $ R 0 
STEP_ 0000 B 0 
STEP_9 00000108 R 0 
TEST_SIGN 00000117 R 0 
TwO_ExP_55 00000000 R 02 
UNDERF LOW 00000122 R 02 
bomen eee wren acanwecce & 
; ;_Psect synopsis : 
PSECT name Allocation PSECT No. Attributes 
. ABS . 00000000 : 0.) 00 ¢ 0.) NOPIC USR CON ABS _ LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
SABSS ess 0.) O1¢ +41.) NOPIC USR CON ABS LCL NOSHR' EXE RD WRT NOVEC BYTE 
_MTHSCODE 000001 311.) 02 ¢ 2.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 
duemenenimmnamMnme 
H Performance indicators ! 
Phase Page faults CPU Time Elapsed Time 
Initialization 34 00:00:00.09 00:00:01.54 
zounend processing 117 00:00:00.45 8: :00:03. 68 
122 00:00:01.45 0:00:05 3 
Synbol table sort 0 00:00: $.93 Se Be 1 
pees ¢ sand Bi % 00:00: 71 00:00:00 a 
m adle outpu H : e : : 
Brect synopsis = 09: 6; 0:03 00:00: 0 § 
Cross-reference output ° 8 :00: 3-9 00:00:00 ° 
Assembler run totals 35 0:00:02.79 00:00:14. 


The working is Limit was 900 pages. 

6542 bytes ( pages) of virtual memory were used to buffer the intermediate code. 

There were 10 pages of symbol table space ot scared to hold 48 non-local and 6 local symbols. 
1 source has were read in Pass 1, producing | object records in Pass 2 

8 pages of virtual memory were used to define 7 macros 


3) 


SSEP=HSES FTiDB:82 EATANTE. SRESArnDmoDeman:1 | P88* (3), 


BTwspnoe 
VAK-11 Macro Run Statistics 


goss es se st es can en ese earn aanae a $ 


Macro Library statistics ! 


' 
énne weer een neem ew me sce nw esse = + 


. 
| 

Macro Library name Macros defined 
-$255$0UA28: CSYSLIBISTARLET.MLB; 2 4 | 
88 GETS were required to define 4 macros. 
There were no errors, warnings or information messages. 
| 

| 


MACRO/ENABLE=SUPPRE SSION/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHDMOD/OBJ=OBJ$:MTHDMOD MSRC$:MTHDMOD/UPDATE=(ENH$:MTHDMOD) 
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